Spin and chiral orderings of frustrated quantum spin chains 
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Ordering of frustrated S—1/2 and 1 XY and Heisenberg spin chains with the competing nearest- and 
next-nearest-neighbor antiferromagnetic couplings is studied by exact diagonalization and density- 
matrix renormalization-group methods. It is found that the S—l XY chain exhibits both gapless 
and gapped 'chiral' phases characterized by the spontaneous breaking of parity, in which the long- 



range order parameter is a chirality, Ki- 



algebraically or exponentially, 
the Heisenberg chains. 
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Such chiral phases are not realized in the S—1/2 XY chain nor in 



Ordering of frustrated quantum spin chains has at- 
tracted considerable interest since these systems exhibit 
a rich variety of magnetic phases due to the interplay be- 
tween quantum effect and frustration. We consider here 
an anisotropic frustrated quantum spin chain described 
by the XXZ Hamiltonian, 



H = 2_^ I jp 2_^ [SiS'g^p -t- s^ Sj. 
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where S^ is the spin-S* spin operator at the £th site, 
Jp > is the antiferromagnetic interaction between the 
nearest- neighbor (p—l) and the next-nearest-neighbor 
(/3=2) spin pairs, and A (0 < A < 1) represents an ex- 
change anisotropy. Note that A = and A = 1 correspond 
to the XY and Heisenberg chains, respectively. 

The ground state phase diagram of the corresponding 
S—1/2 system has been extensively studied either nu- 
merically [P|] or analytically l^^^l- These studies have 
revealed that when J2 is smaller than a critical value, 
i.e., j= J2/ Ji < ic the system is in the gapless spin-fluid 
phase in which the antiferromagnetic spin correlation de- 
cays algebraically. By contrast, for larger values of j > jc, 
the system is in the dimer phase with a finite energy gap 
above the doubly degenerate ground states. The dimer 
phase is characterized by the spontaneously breaking of 
both parity and translation symmetries with preserving 
time-reversal symmetry. The value of jc has been esti- 
mated to be jc ^.241 for the Heisenberg chain y]. Al- 
though there is no magnetic long-range order (LRO), the 
nature of the magnetic short-range order (SRO) changes 
at the Lifshitz point j^, where ji~0.5 for the Heisen- 
berg chain W. For j < jl, the system has the standard 
Neel-type antiferromagnetic SRO and the structure fac- 
tor S{q) has a maximum at q — n, while for j > ji^, the 
system has a helical SRO with the maximum of S{q) at 
some q=Q < tt. 

In the case of S—1, by contrast, no dimer phase oc- 
curs P,B| . The Heisenberg chain is in the Haldane phase 
characterized by a singlet ground state and a finite en- 



ergy gap above it. A first-order transition takes place at 
j—JT — 0.744 between the 'single-chain' Haldane phase 
£^t j < JT and the 'double-chain' Haldane phase at j > jt 
Q. In the XY case, on the other hand, the situation re- 
mains not entirely clear. Analytical studies based on the 
bosonization method suggested that the gapless phase at 
j = (the so-called XYl phase) extended to finite j > 
whereas numerical studies suggested that the Haldane 
phase was stabilized for J > [p[. In any case, the fate 
of such XYl 01 Haldane phase at larger j has not been 
clarified. 

In the classical limit S* — > 00, the system exhibits 
a magnetic LRO, either of the Neel-type {j < 1/4) 
or of the helical-type characterized by the wavenumber 
g = cos"^(-l/4j) (j > 1/4). In the XY case, such heli- 
cally ordered state possesses a twofold discrete degener- 
acy according as the helix is either right- or left-handed, 
in addition to a continuous degeneracy associated with 
the original U{1) symmetry of the XY spin. This dis- 
crete degeneracy is characterized by mutually opposite 
signs of the total chirality defined by M 
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Si Sf_^_i — [Sj X Si+i] 



i+l\z, 



where N is the total number of spins. Chirality is in- 
variant under both U{1) spin- rotation and time-reversal 
operations, but changes its sign under parity operation. 
In the Heisenberg case, while there no longer exists a 
discrete chiral degeneracy, one can still formally define 
the chirality by Eq.(p|) as a 2-component of the vector 
chirality, Ki=SjxSi-|-i. Note that the chirality defined 
above is distinct from the scalar chirality of the Heisen- 
berg spin often discussed in the literature pQI defined 
by Xi — Si_i • Sj X Si-|_i: The scalar chirality of the 
Heisenberg spin changes sign under the time-reversal op- 
eration unlike the chirality (2). Since the classical chain 
(1) always has a planar spin order, the scalar chirality Xi 
vanishes trivially. 



Recent studies on various frustrated classical systems 
with continuous symmetry have revealed that such chiral 
degrees of freedom often give rise to novel ordering be- 
haviors such as phase transitions of new universality class 
pj| , novel magnetic phase diagrams with new multicriti- 
cal behavior [0, and even a novel 'chiral phase' in which 
only the chirality exhibits a LRO without the standard 
spin LRO p2|-p^. Meanwhile, systematic studies of the 
possible chiral order in frustrated quantum spin chains 
has been scarce so far. 

In this Letter, we perform a numerical study of the 
ground-state properties of a class of S—1/2 and 1 frus- 
trated spin chains (1) in search for a possible chiral or- 
der. We have found that both gapless and gapped chiral 
phases, where the chirality has a finite LRO while the 
spin correlation falls either algebraically or exponentially, 
are realized for a wide range of j in the S=l XY chain, 
but not in the S=l/2 XY chain nor in the 5=1/2 and 1 
Heisenberg chains. 

In order to probe the possible chiral order, we first 
calculate the square of the total chirality < k^ >, where 
the chiral order parameter is defined by Eq.(2) now for 
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(b) S=1 XY 
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FIG. 1. Binder parameter of the chirality Qk. versus j: (a) 
5=1/2 Xychain; (b) S=l J^ychain. 



quantum spins, together with the associated Binder pa- 
rameter, 
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(3) 



by exactly diagonalizing finite open chains with even N 
up to iV=20 (S'=l/2) or 7V=16 (5=1) (l|. All ground 
states are found to belong to the subspace of S^^^^y = 
with even parity. 

The calculated Binder parameters of the chirality 5^ 
of the XY chain are shown in Fig.l for various A^ as a 
function of j. As shown in Fig. 1(a), g^, of the 5 = 1/2 
chain constantly decreases with increasing A^ for any j, 
demonstrating the absence of chiral order. In contrast, 
(7k of the 5=1 chain exhibits a significantly different 
behavior. As can be seen from Fig. 1(b), g,^ for various N 
cross at almost the same j, j = 0.475±0.005, above which 
(7k increases with increasing N, indicating the existence 
of a finite chiral LRO for larger j. 

In the case of Heisenberg chains, on the other hand, 
we have found that for both cases of 5 = 1/2 and 1 the 
calculated g^ (not shown here) constantly decreases with 
increasing A^, which clearly shows that the chiral order 
is not realized in the Heisenberg chains. 

In order to examine more closely the nature of the 
transition and of the possible phases, we calculate, con- 
centrating on the 5=1 XY open chain, the two-point 
chiral, spin, and string correlation functions defined by 



*-'k('") — < '*ro-r/2'^ro+r-/2 >7 

ro+r/2-1 
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j=ro-r/2+l 
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^j)^ro+r/2 >' (6) 



by means of the density-matrix renormalization-group 
method. [Em. Here ro represents the center position of 
open chain, i.e., ro—N/2 for even r and ro=(A^-l- 1)/2 for 
odd r. We employ the infinite system method by using 
Af block states {M < 300) in the subspace of S^^^^y=0 
with even parity. Convergence of the results with respect 
to M has been checked by consecutively increasing M. 
In the chiral phase, the chiral correlation function Cuir) 
tends to a finite constant at large r, while in the Haldane 
phase, the string correlation function Cstr{T) tends to a 
finite constant. 

The calculated r-dependence of the chiral, string, and 
spin correlation functions are shown in Fig.2(a)-(c) on 
log- log plots for several typical values of j. As can be 
seen from Fig. 2(a), the data of C^{r) for j > jd — 0.473 
bend up at larger r suggesting a finite chiral LRO, while 
they bend down for j < jd suggesting the absence of chi- 
ral order. Our estimate of jd — 0.473 ± 0.001 is consis- 
tent with the estimate from the Binder parameter given 
above. Meanwhile, as shown in Fig. 2(b), the data of 



Cstrir) for j < jc2 — 0.490 bend up for larger r sug- 
gesting a finite string order characteristic of the Haldane 
state, while for j > jc2 they show a linear behavior sig- 
naling a power-law decay of the string correlation. The 
existence of the Haldane phase for smaller j is consistent 
with the previous finding of Ref. M . 

An interesting observation here is that our estimate 
of jc2 = 0.490^gQ|!)5 is distinctly larger than that of 
jci = 0.473 ± 0.001, which means that there exist two 
different types of chiral phases ^ one with a finite string 
order (jd < j < jc2) and the other without the string 
order (j > jc2)- This finding is also supported by the be- 
havior of the spin correlation function Cs{r) in Fig. 2(c), 
in which Cs{r) divided by the oscillating factor cos{Qr) 
is shown. Indeed, the data of Cs{r) exhibit a linear be- 
havior for j > jc2 indicating a power-law decay of he- 
lical spin correlations (gapless state), while for j < jc2 
they bend down suggesting an exponential-decay of heli- 
cal spin correlations (gapped state). Note that above the 
Lifshitz point, which we estimate to be jl = 0.313±0.001 
for the 5=1 XY chain, the system exhibits a helical 
SRO characterized by a wavevector Q < tt. On the other 
hand, the absence of magnetic (spin) LRO has rigorously 
been proven for any j and for general S < oo ||l^ . The 
existence of a novel intermediate phase, a gapped chiral 
phase, in a narrow but finite range of j can clearly be seen 
in Fig. 2 from the behavior of the correlation functions at 
j — 0.477 which lies between jci and jc2- 

Thus, on increasing j, the S = 1 XY spin chain 
undergoes two successive transitions, first at j = jd 
from the Haldane phase with no chiral order to the 
gapped chiral phase (or the chiral Haldane phase), then 
at j = jc2 from the gapped chiral phase to the gapless 
chiral phase. In the gapped chiral phase, the chiral and 
string LRO coexist with exponentially-decaying spin cor- 
relations, whereas in the gapless chiral phase, only the 
chirality shows a LRO with algebraically-decaying spin 
and string correlations. In the gapped (gapless) chiral 
phase, the ground state is doubly degenerate, each of 
which is characterized by the opposite sense of the chi- 
rality, with (without) a finite gap above it. 

The gapped chiral phase is characterized by a spon- 
taneously broken parity with preserving both transla- 
tion and time-reversal symmetries. From a broken sym- 
metry, the transition at j = jd is expected to be of 
the Ising-type. We extract the chiral correlation length 
^K from the calculated CK{r) and fit it to the standard 
power-law form ^^ '^ (jci ~ j)'^" ■ Our present estimate, 
i/k — 0.9 ± 0.1 appears to be slightly smaller than, but 
not inconsistent with the 2D Ising value v^ — 1. Further 
detailed study of the critical properties, including the na- 
ture of the KT-like transition at j = _jc2, is in progress. 

It should be noticed that Nersesyan et. al. recently 
discussed for the S'=l/2 XY chain in the limit of large j 
the possibility of the parity breaking jl8| . In contrast to 
this suggestion, for the S=l/2 XY and Heisenberg 
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(c) spin correlation 
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FIG. 2. Correlation functions versus r on log-log plots for 
various j: (a) chiral correlation CK{r); (b) string correlation 
— Cstr(r); (c) spin correlation Cs{r) divided by the oscillating 
factor cos{Qr). Arrows in the figures represent the extrapo- 
lated r — oo values. The number of block states is equal to 
M = 300. To illustrate the M-dependence, we also indicate 
by crosses the data for M — 220 and 260 for several cases 
where the M-dependence is relative large. Note that some of 
the data points for larger r are omitted for clarity. 

chains, we did not find the chiral phase p9| | nor any 
new phase except for the well-known spin-fluid and dimer 
phases, corroborating the type of the phase diagram pre- 
viously reported by Haldane H] and by Tonegawa et al 



y. The chiral phase found here appears to be specific 
to the 5' = 1 XY chain. Such tendency may roughly be 
understood by noting that only the XY chain sustains an 
Ising-like discrete chirality which has a stronger ordering 
tendency than the continuous spin variable while quan- 
tum fluctuations might be strong enough in the S = 1/2 
chain to wash out even the chiral ordering. 

It might also be interesting to notice that the observed 
novel transition behavior of the S = 1 XY chain has a 
similarity to that of the 2D frustrated classical XY mod- 
els such as the triangular-lattice XY antiferromagnet p^ ] 
or the Josephson-junction array in a magnetic field p3|. 
In these classical systems, Miyashita and Shiba, and more 
recently Olsson observed by Monte Carlo simulations 
that the thermal phase transitions occurred in two steps 
with two types of chiral phases, each characterized by 
exponentially-decaying and algebraically-decaying spin 
correlations. At the moment, we donot know whether 
there exists a deeper connection between the two sys- 
tems. 

Finally, we wish to briefly discuss the possible experi- 
mental implication of our results. In order to observe the 
chiral phase, one needs to prepare an 5* = 1 XY zig-zag 
chain with its j value in a suitable range. In the pres- 
ence of weak 3D interchain interaction, while the gapless 
chiral phase stabilized for larger j is expected to show 
the standard helical spin LRO, the gapped chiral phase 
could remain gapped. Hence, it is challenging to experi- 
mentally observe the gapped chiral state in an appropri- 
ate model material. One problem here might be that the 
gapped chiral phase is realized in a rather narrow range 
of j. It might then be necessary to tune the j value by 
some experimental method, such as by applying pressure. 
Once an appropriate sample could be prepared, it is in 
principle possible to measure the chirality by using, e.g., 
polarized neutrons. [ [ll|J2(]|] 

In summary, from numerical studies of the ground- 
state properties of a class of 5= 1/2 and 1 frustrated spin 
chains we have found that both gapless and gapped chi- 
ral phases, where the chirality has a finite LRO while the 
spin correlation falls either algebraically or exponentially, 
are realized for a wide range of j in the S'=l XY chain, 
but not in the S'=l/2 XY chain nor in the 5=1/2 and 1 
Heisenberg chains. Further details of the results includ- 
ing the critical properties and the full phase diagrams 
will be reported elsewhere. 

We thank Prof. T. Tonegawa for valuable discussion. 
Numerical calculations were carried out in part at the 
Yukawa Institute Computer Facility, Kyoto University. 
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